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Abstract 

We consider the implications from the possibility that the recently observed state X(3872) is a meson- 



antimeson molecule. We write an effective Lagrangian consistent with the heavy-quark and chiral sym- 

o _ 

metries needed to describe X(3872). We claim that if X(3872) is a molecular bound state of D*° and D° 
r~| ■ mesons, the heavy-quark symmetry requires the existence of the molecular bound state Xf, of B*° and B° 

a: 

' , with the mass of 10604 MeV. 
Ph. 



PACS numbers: 12.39.Hg, 12.39.Fe, 12.39.Mk, 14.40.Gx 
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I. INTRODUCTION 



The past two years have seen several experimental observations of new heavy and light quark 
states. Most of those states contain a charm quark, so their observation rekindled interest in 
heavy- flavor spectroscopy The unusual properties of those states invited some speculations 
regarding their possible non-gg nature. Among those is the X(3872) state which, being discovered 
in the decay X(3872) — > J/ip-n + -n~ , contains charm- ant icharm quarks . While a traditional 
cc quarkonium interpretation of this state has been proposed m, its somewhat unusual mass and 
decay patterns prompted a series of more exotic interpretations p. Since the mass of the X(3872) 
state lies tantalizingly close to the D*°D° threshold of 3871.3 MeV, it is tempting to assume that 
X(3872) could be a D*°D° molecular state 0,0. Recent Belle data appear to be consistent with 
this assignment, preliminarily confirming its J PC = 1 ++ quantum numbers js). Of course, states 
of different "nature" can mix, if they have the same quantum numbers, further complicating the 
interpretation of the experimental data p 1 . 

An unambiguous identification of this state must be done with many different measurements of 
its decay and production patterns. Regardless of whether A (3872) is identified to be a molecule 
or a regular qq charmonium, a theoretical analysis of heavy-meson molecular states should be 
performed. Until recently |7[ these studies were done mostly with the help of various quark 
models 0, Q|. In this paper we shall study those states using the techniques of effective field 
theories. 

This study is possible due to the multitude of scales present in QCD. The extreme smallness of 
the binding energy 

E h = (m D0 +m D0 ,) - M x = -0.6 ± 1.1 MeV (1) 

suggests that this state can play the role of the "deuteron" (or "deuson," see N. A. Tornqvist's 
paper in (|) in meson- meson interactions. The "deuteron-like" nature of this state allows us to use 
methods similar to those developed for the description of the deuteron, with the added benefit of 
heavy-quark symmetry. The tiny binding energy of this molecular state introduces an energy scale 
which is much smaller than the mass of the lightest particle, the pion, whose exchange can provide 
binding. Thus, for a suitable description of this state in the framework of effective field theory, 
the pion, along with other particles providing possible binding contributions (i.e. the p-meson and 
other higher mass resonances), must be integrated out. The resulting Lagrangian should contain 
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only heavy-meson degrees of freedom with interactions approximated by local four-boson terms 
constrained only by the symmetries of the theory. This approach is similar to the description of 
the deuteron in the effective theory without dynamical pions 0] . Nevertheless, we shall often 
appeal to the "exchange picture" to gain insight into the structure of the effective Lagrangian. 

This approach provides a model-independent description of molecular states with somewhat 
limited predictive power. In particular, we would not be able to say whether the state V(3872) 
is a D*°D° molecule or not. What we would be able to say is that if indeed X(3872) is a D*°D° 
molecule, the heavy-quark symmetry makes a definite statement on the existence of a molecular 
state in the B*°B° system. We also show that even though D and D* are degenerate in the 
heavy-quark limit, the existence of a molecular state in D*°D° channel does not necessarily imply 
a molecular state in the D°D° or B*°D° channels. 

This paper is organized as follows. In Sec. II we write the most general effective Lagrangian 
consistent with the heavy-quark and chiral symmetry. In Sec. Ill we obtain the bound-state 
energy by solving a system of Lippmann-Schwinger equations and relate the bound-state energies 
of D*°D° and B*°B° states. We present our conclusions in Sec. IV. 



II. THE EFFECTIVE LAGRANGIAN 



In order to describe the molecular states of heavy mesons we need an effective Lagrangian which 
contains two- and four-body interaction terms. The two-body effective Lagrangian that describes 
the strong interactions of the heavy mesons P and P* (P = B, D) containing one heavy quark Q 
is well known jl^ : 

1 



C, 



+ ^Tr 



H {Q) v ■ DHM 



2m i 



-Tr 



fJbV 



- £<9^ 1 



+ ... 



(2) 



where the ellipsis denotes terms with more derivatives or including explicit factors of light quark 
masses, D^ b = <5 a 6<9 M — (1/2) (^<9 M £ + £<9 M £"n , and g is the P*Pn coupling. As usual, we introduce 

V 'ah 

a superfield describing the combined doublet of pseudoscalar heavy-meson fields P a = (P°, P + ) 
and their vector counterparts with v ■ P£ = 0, 

1+ jb 



-jj(Q)a 
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°pfi Q )y. 



(3) 



These fields have the usual transformation properties under the heavy-quark spin symmetry and 
SU(2)y flavor symmetry 1 , 

H (Q) ^ S (HMrf) a , H (Q)a - (UH {Q) J S~\ (4) 

and describe heavy mesons with a definite velocity v (l4|. The third term in Eq. (J2J) is needed to 
account for the P — P* mass difference A = m P » — m P = —2\2/m P . 

The pseudo-Goldstone fields are introduced as £ = e lM ^, where M is the usual meson matrix 



M 



(5) 



and / ~ 135 MeV is the pion decay constant. Notice that since heavy quark- ant iquark pair 
production is absent in this effective theory, the effective Lagrangian of Eq. (j2J) does not contain 
heavy antimeson degrees of freedom. Since we are describing the molecular states of heavy mesons, 
those fields should have to be explicitly added to the Lagrangian. The fields Hfp> and Hjp" 
that describe the propagation of heavy antimesons, i.e. containing the heavy antiquark Q, are 
introduced as 

HP = [P;f V - Pp 75 ] = l°Hph°, (6) 

and transform as HP -> (UH®)) a S- 1 and H iQ)a -> S (H {Q) Wy under heavy-quark spin and 
SU(2)y symmetry. 

In order to write an effective Lagrangian describing the properties of X(3872), we need to couple 
the fields HjQ) and H^ a so that the resulting Lagrangian respects the heavy-quark spin and chiral 
symmetries. Since the binding energy of X(3872) is small, the size of a bound state is rather large. 
This means that the particular details of the interaction of the heavy meson and antimeson pair 
(for example, a p-meson exchange contribution) are irrelevant for the description of such a bound 
state and can be well approximated by four-meson local interactions. One can write a Lagrangian 
describing X(3872) by first writing an effective Lagrangian above ji = and then matching it 
onto the Lagrangian for fi < m n , i.e. integrating out the pion degrees of freedom. 

The general effective Lagrangian consistent with heavy-quark spin and chiral symmetries can 
be written as 

C = C 2 + £ 4 , (7) 



1 A generalization of this discussion to the flavor SU(3) symmetry is rather straightforward. 
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where the two-body piece, consistent with reparametrization invariance, is given by Eq. (J2J) and 
the four-body piece is 



£ 4 =^Tr 



fj(Q) H (Q) 



7n 



Tr 



r 



+ -±Tr 

4 



7^75 



Tr 



7 M 75 



This Lagrangian, together with £2, describes the scattering of P and P* mesons at the energy scale 
above m^. Integrating out the pion degrees of freedom at tree level corresponds to a modification 
C 2 — > C 2 + (2/3) (g/ f) 2 . Since in this paper we will not discuss matching at higher orders, the 
Lagrangian in Eq. (jEJ) will be used for the calculation of the bound state properties of X(3872). 

By virtue of the heavy-quark spin symmetry, the same Lagrangian governs the four-boson 
interactions of all P(*> = D^*' or B^ states, while the flavor SU(2)v implies that there could 
be four such states for each P^P^\ Indeed, not all of these states are bound. Here we shall 
concentrate on X(3872), which is a bound state of two neutral bosons, P a = P° = P, assuming 
the isospin breaking advocated in Ref. p. Notice that the most general Lagrangian involves two 
couplings, C\ and C<i- Other Dirac structures are possible, but will yield the same Lagrangian for 
the PP* bound state. 

In order to relate the properties of PP* molecules in the charm and beauty sectors we shall need 
to see how C{ couplings scale as functions of the heavy-quark mass M. To see this, we recall that a 
system of two heavy particles requires a nonrelativistic v/c expansion, not a 1/M expansion. This 
is necessary to avoid that the resulting loop integrals acquire pinch singularities [llj. Therefore, 
we must powercount p° ~ jf/M, where p is a characteristic 3-momentum of a heavy meson in the 
PP* molecular bound state, which implies that the first and the second terms in Eq. (j2J) scale in 
the same way. Since the action S = J d A x L does not scale with the heavy-quark mass, this implies 
that d 4 x ~ M and the Lagrangian density £ ~ 1/M. The kinetic term in Eq. (|2j) then gives the 
expected scaling of the heavy- meson field H ~ P ~ P* ~ M°, which in turn implies from Eq. (jSJ) 
that the couplings 

a ~ i/M. (9) 

This dimensional analysis, however, cannot be used to predict the relative contributions of other 



couplings in £ 4 , say, relativistic corrections to Eq. 
molecular state close to threshold in the first place 



because of the fine-tuning which produces a 



We will use it only to relate properties of 



DD* and BB* systems. Similar dimensional analysis was proposed for non-relativistic QCD 



12]. 
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Evaluating the traces yields for the PP* sector 

£ 4 pp , =-C 1 p(Q)t p(Q) p*(Q)t p*(Q)v _ (7 1 p*(Q)t p*(Q)A< p(Q)t p(Q) 

+C 2 P (Q)t P*( Q )p*^ )tM P® + C 2 P* (Q)t p( Q) p(<2)t p*(Q)M + . . . (10) 

Notice that this Lagrangian differs from the one used in jjj, where the interaction strength is de- 
scribed by a single parameter A = —C\ = C 2 . The difference can be understood in the "exchange" 
model, where C\ and G% come from the exchanges of mesons of different masses and parity. In 
this language the model of Ref. j?| corresponds to the situation of degenerate parity states. In 
QCD, however, negative parity states are generally lighter than their positive parity counterparts. 
This is especially clear for the lightest octet of pseudoscalar mesons, where chiral symmetry forces 
their masses to be almost zero, while all the corresponding scalar mesons have masses of the order 
of 1 GeV. We will nevertheless show that the resulting binding energy still depends on a single 
parameter, a linear combination of C\ and C 2 . 

Similarly, one obtains the component Lagrangian governing the interactions of P and P, 

£ 4PP = c 1 pW)tp(«)p(Q)tp(Q). (11) 

Clearly, one cannot relate the existence of the bound state in the PP* and PP channels, as the 
properties of the latter will depend on G\ alone, not a linear combination of C\ and C 2 . 



III. PROPERTIES OF BOUND STATES 

In order to describe bound states we shall modify the approach of S. Weinberg The lowest- 
energy bound state of P and P* is an eigenstate of charge conjugation. Here we carry out our 
calculation for (P,P*) = (D,D*), but analogous considerations will apply to the B system. The 
two eigenstates of charge conjugation will be given by 

1 



l*±> 



V2 



\D*D) ± \DD ) 



(12) 



To find the bound-state energy of X(3872) with J PC = 1 ++ , we shall look for a pole of the 
transition amplitude T ++ = (X + \T\X + ). 

We first define the following transition amplitudes, 



T n ={D*D\T\D*D), T 12 = (D*D\T\DD*), 




D* 



D 











FIG. 1: Transition amplitudes for the D — D scattering written in the form of Lippmann-Schwinger 
equations. Double lines indicate the vector D* or D* mesons, solid lines pseudoscalar D or D states. 

T 21 = (DD*\T\D*D), T 22 = (DD*\T\DD*), (13) 

which correspond to the scattering of D and D* mesons. Clearly, at tree level, Ta ~ C% and 
Tij ~ C 2 , since we consider only contact interactions. But we also have to include a resumma- 
tion of loop contributions to complete the leading order [ll|. These transition amplitudes satisfy 
a system of Lippmann-Schwinger equations depicted in Fig. [TJ 

iT lx =-iC x + J -jL*pT xl Gpp.C x - J -^Ti2G PP *C 2 , 
d^q f d A q 



iT u = iC 2 -fj^T n Gpp.C a + f^T 3a Gpp.C 1 , 
iT 21 = %C 2 + J A^T^Gpp.d - J -0^T 22 G PP *C 2 , (14) 



iT, 



■>■-> - — zCi — / —ToiGpp*Co ~\~ / ~, —TooG pp*Ct , 

J (2tt) 4 J (2vr) 4 



where 



Gpp* — — r— r-, (15) 

4 (p 2 /2m D , + q - A - q 2 /2m D , + iej [p 2 /2m D — q — (p/2m D + iej 
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p is the momentum of one of the mesons in the center-of-mass system, and we canceled out factors 
of m D m D , appearing on both sides of Eq. (|14jl. Note that in Eq. (JT5)l the vector propagator 
includes the mass difference A, as a consequence of the term proportional to A2 in the Lagrangian 
of Eq. (0) [13]. This choice of propagators is not unique, but our results will not depend on 
it, because it amounts to a choice of a finite phase multiplying the heavy- meson fields. This 
rephasing is equivalent to measuring energies with respect to the pseudoscalar mass, mp. Then, 
the position of the transition amplitude pole on the energy scale will be measured with respect to 
the "constituent mass" of the system, which is in our case twice the pseudoscalar mass m D \\\ - 
A different choice of phase will give different propagators but also a different "constituent mass" . 

Since we are interested in the pole of the amplitude T ++ , we must diagonalize this system of 
equations rewritten in an algebraic matrix form, 



1 11 



T 12 
\T22J 



C 2 



( 



+ iA 



\ 



V 



(rr \ 



(16) 



-Ci c 2 
C 2 -d 
-d c 2 
C 2 -dj 

Notice that the matrix is in the block-diagonal form, which allows us to solve Eq. (fTB*j) in two steps 
working only with 2x2 matrices. The solution of Eq. |T6j) produces the T ++ amplitude, 



T 21 
\T 22 ) 
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Tl2 + Tn\ + T< 



A 



-22 



iXA 



(17) 



with A = —C\ + C 2 , and A is a (divergent) integral 



A 



1/ 



(2vr) 4 (p>/2m D , + q - A - q*/2m D « + ie) (f/2 
d 3 q 1 



D 



q - cp/2m D + ie 



(27r) 3 q 2 - 2[i DD « (E — A) — ie " ' 

where E = p*/2fi DD ,, fi DD * is the reduced mass of the DD* system, and we have used the residue 
theorem to evaluate the integral over q°. The divergence of the integral of Eq. (|18|) . as usual, is 
removed by renormalization. We choose to define a renormalized A^. within the MS subtraction 
scheme in dimensional regularization. In this scheme the integral A is finite, which corresponds 
to an implicit subtraction of power divergences in Eq. (|18|). Computing the second integral in 
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Eq. (fTSjl by analytically continuing to d — 1 dimensions yields 

A, 

(19) 



A = -tt^ddApI 



This implies for the transition amplitude 

A 



r3 



T ++ = " R (20) 

1 + (z/8vr)A R \p\Jl-2fi DD ,A/jP 



The position of the pole of the molecular state on the energy scale can be read off Eq. (j2U|l . 



32tt 2 

^R^DD* 



Ep l e = ~2 — 3 A. (21) 



This is the amount of energy we must subtract from the "constituent mass" of the system, deter- 



mined above as 2m D , in order to calculate the mass 



M x = 2m D - E Pole = 2m D + A - — . (22) 



32tt 2 

^rIAdD* 

Recalling the definition of binding energy, Eq. and that m D * = m D + A, we infer 

^R^DD* 

Assuming Ef, = 0.5 MeV, which is one sigma below the central value in Eq. and the 

experimental values for the masses 



17| . we obtain 

\ R ~ 8.4 x 10~ 4 MeV^ 2 . (24) 

Note that the binding energy scales as 1/M in the heavy-quark limit. Thus, the smallness of the 
binding energy is implied in the heavy-quark limit. The small binding energy of the X(3872) state 
implies that the scattering length a D is large and can be written as 



(2^ DD ,E b y 1 = (25) 



yielding a numerical value a D = 6.3 fm. Since the scattering length is large, universality implies 
that the leading-order wave function of X(3872) is known, 

e -r/a D 

4>DD*{r) = ^— (26) 
27ia D r 



This can be used to predict the production and decay properties of X(3872) tIQ]. 

Once we establish the molecular nature of X(3872), its experimental mass gives us its binding 
energy. The latter is dependent on the coupling constant \ R , and may be used to predict the 
binding energies of hypothetical molecular bound states in the beauty sector, as well as to discuss 
its implications for the beauty-charm sector. Alternatively, the coupling constant A^ can be fixed 
from the resonance-exchange model, in which case (model-dependent) predictions are possible for 
all the heavy sectors. Taking into account the scaling of A^ with M given by Eq. ©, we obtain 

X R ~ X R • i 2 ') 

r'BB* 

Formula (|2*3*j) can now be used for the B system. This implies the existence of an S-wave bound 
states with binding energy Ef, = 0.18 MeV, mass Mx b = 10604 MeV, and a scattering length a B 
= a D = 6.3 fm, because of the lack of scaling of the scattering length with the heavy-quark mass. 
The above prediction for the binding energy is lower than the quark-model prediction of Wong 
in Similar considerations apply to D°D° and B°B° states: In their case the starting point is 
the Lagrangian term in Eq. (jllj) . Since it involves only a single term, the calculations are actually 
easier and involve only one Lippmann-Schwinger equation. The resulting binding energy is then 

F 256?r2 (?R) 
L b - r2 3 , (28) 

and the equivalent formula for the B°B° system may be obtained by rescaling the coupling constant 
Cf R as we did with X R in Eq. (|27|). Examining Eq. (|2*Hj) we immediately notice that the existence 
of a bound state in the D*D channel does not dictate the properties of a possible bound state in 
the D°D° or B°B° channels, since C% and C2 are generally not related to each other. 

The discussion of the beauty-charm system parallels what was done above. In this case the 
situation is a bit different, because there are two states, B°D*° and D°B*° . The treatment of 
these states depends on how the heavy-quark limit is taken. They have different masses, since 
Abb* 7^ Add*. This implies that these two states do not mix and have to be treated separately, so 
that their binding energies could be obtained by respectively substituting A^ — > C\r and fi D D* ~^ 
fi BD * or Hdb* i n Eq. (|23jl . But just like in the case of D°D° or B°B° channels, we cannot predict 
bound states in these channels from the heavy-quark symmetry arguments alone. 
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IV. CONCLUSION 



We have used an effective field theory approach in the analysis of the likely molecular state 
X(3872), by describing its binding interaction with contact terms in a heavy-quark symmetric 
Lagrangian. The flexibility of this description allows us to ignore the details of the interaction and 
to concentrate on its effects, namely a shallow bound state and a large scattering length. Taking 
into account the universality and the scaling of the effective coupling constants we are able to 
extend our description to the B system. We found that if X(3872) is indeed a molecular bound 
state of D*° and D° mesons, the heavy-quark symmetry requires the existence of the molecular 
bound state X b of B*° and B° with the mass of 10604 MeV. 
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